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We study the wetting behaviour of a symmetrical binary fluid below the demixing temperature 
at a non-selective attractive wall. Although it demixes in the bulk, a sufficiently thin liquid film 
remains mixed. On approaching liquid/vapour coexistence, however, the thickness of the liquid 
film increases and it may demix and then wet the substrate. We show that the wetting properties 
are determined by an interplay of the two length scales related to the density and the composition 
fiuctuations. The problem is analysed within the framework of a generic two component Ginzburg- 
Landau functional (appropriate for systems with short- ranged interactions). This functional is 
minimized both numerically and analytically within a piecewise parabolic potential approximation. 
A number of novel surface transitions are found, including first order demixing and prewetting, 
continuous demixing, a tricritical point connecting the two regimes, or a critical end point beyond 
which the prewetting line separates a strongly and a weakly demixed film. Our results are supported 
by detailed Monte Carlo simulations of a symmetrical binary Lennard- Jones fluid at an attractive 
wall. 

PACS numbers 68.45.Gd, 68.10.-m, 05.70.Jk, 68.15.-Fe 
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Phase transitions are well known to be influenced by 
geometrical confinement In practice, confinement is 
often imposed by rigid external constraints, for example 
the surfaces of porous or artificially nanostructured me- 
dia. However, it can also be an inherent feature of a 
system, as occurs for a liquid wetting film bound to a 
solid substrate and in equilibrium with its vapour |^ . In 
such a situation the liquid is confined between the rigid 
substrate and the flexible liquid-vapour interface. 

The effects of confinement are particularly pronounced 
in the region of critical points. Under such conditions the 
system exhibits strong order parameter fiuctuations, the 
correlation length of which may become comparable with 
the linear dimension of the confined system. When this 
occurs the effects of confinement are felt not just near the 
confining surfaces, but propagate throughout the system 

Critical fiuctuation are relevant to the properties of liq- 
uid wetting films if the liquid in question possesses an ad- 
ditional internal degree of freedom. Then the state of the 
liquid is described not only by its number density, but by 
an additional parameter measuring the degree of internal 
order. Examples are binary liquids, where the additional 
order parameter is the relative concentration of species, 
and ferrofiuids where it is the magnetisation. In such 
systems the geometrical constraint (i.e. the film thick- 



ness) can itself depend on the state of order in the liquid 
film. For example, simulation and experiment have re- 
cently shown that critical concentration fiuctuations can 
change the equilibrium thickness of a wetting layer of a 
binary liquid — the so-called critical Casimir effect 1^-^ . 

It transpires, however, that the interesting conse- 
quences of interplay between the degree of order in the 
wetting layer and film thickness are not limited to the 
critical point itself. To illustrate this, it is instructive 
to consider the following Gedanken experiment. Let us 
take a symmetrical binary fluid, i.e. a fluid in which 
particles of the same species have one strength of in- 
teraction, while interactions between dissimilar species 
have another strength. As elucidated in ref. 0, it is 
possible to arrange for such a system to exhibit a line of 
continuous demixing transitions, terminating in a critical 
end point on the liquid side of liquid- vapour coexistence. 
Suppose now that the fluid is placed in contact with a 
non-selective attractive substrate (wall) acting equally 
on both species. If the wall is sufficiently attractive, 
complete wetting occurs at and above the critical end 
point temperature T^ep as liquid-vapour coexistence is 
approached from the vapour side. But what happens for 
T < Tcep"! Far from coexistence, the wetting film is suf- 
ficiently thin that demixing will certainly be suppressed. 
On approaching the coexistence curve, however, the film 



thickness grows and it is tempting to argue that it even- 
tuaUy exceeds the correlation length of composition fluc- 
tuations, whereupon the film spontaneously demixes. 

Notwithstanding the appealing simplicity of this argu- 
ment, it turns out to contain two flaws which render the 
actual situation rather more complex. First, the thick- 
ness of the mixed wetting film will not increase beyond all 
limits below the critical end point Tcep- This is because 
a hypothetical mixed bulk liquid would not coexist with 
the vapour phase at the same chemical potential fio as a 
demixed liquid, but rather at a chemical potential which 
is shifted by 



(1.1) 



towards the liquid side in the fJ. — T plane . Since the 
thickness of the mixed film grows as 
is bounded from above by 



/* oc ln(/i* - ^o) « ln(T'cep - T). 



(1.2) 



The maximum thickness l^: of a hypothetical mixed film 
thus diverges logarithmically on approaching Tcep- In 
contrast, the correlation length ^ diverges much faster, 
like ^ oc {Tcep — T)"" and will thus always exceed l^, 
sufficiently close to Tcep- 

The second flaw in our argument is its implicit assump- 
tion that the composition or order parameter profile is 
confined in an effectively steplike density profile, i.e. that 
the interfacial width between the liquid and the vapour is 
much smaller than the correlation length of composition 
fiuctuations. Although this is true in the region of the 
critical end point, for temperature sufficiently below Tcep 
the correlation lengths of density and composition fluc- 
tuations can be comparable and the interplay between 
the two subtle. 

In the present study, we deploy mean fleld calculation 
and Monte Carlo simulation to elucidate the range of 
possible wetting behaviour of a symmetrical binary fluid 
mixture at a non-selective attractive wall for tempera- 
tures below Tcep- For the sake of simplicity, we have 
chosen to ignore long range dispersion forces in the an- 
alytical calculations, instead taking the interactions to 
be short ranged. This allows us to base our study on a 
generic Ginzburg-Landau model which we solve numer- 
ically and analytically within a square gradient approx- 
imation. The latter leads to the construction of a fllm 
free energy (effective interface potential) highlighting the 
role of the different length scales involved in the problem. 
We show that the competition of length scales results in 
wetting phase behavior considerably more complex than 
has hitherto been appreciated. The analytical results 
are compared with (and supported by) detailed Monte 
Carlo simulations of a binary Lennard- Jones fluid in a 
semi-infinite geometry, interacting with a non-selective 
attractive substrate via dispersion forces. 

With regard to previous related work, the sole discus- 
sion of wetting of symmetrical binary fluids at a non- 
selective wall (of which we are aware) is that of Dietrich 



and Schick ||Tl|] who considered them in a sharp kink 
approximation treatment of binary fluids having long- 
ranged interactions. Most other work on the wetting 
properties of binary fluids has focused on the case of a se- 
lective substrate (favouring one component) [p^-p^. Al- 
though such models correspond more closely than ours to 
experimental conditions jl^Jl^, they lack the aspect of 
simultaneous demixing/ordering and wetting which is of 
interest to us here. It should be stressed, however, that 
realisations of fluids having symmetrical internal degrees 
of freedom do in fact exist, notably in the form of ferroflu- 
ids , so our model is of more than purely theoretical 
interest. 

More general studies of wetting in systems with more 
than one order parameter and associated length scales 
have been discussed by Hauge ||l^ , who pointed out that 
wetting exponents may become nonuniversal even on the 
mean field level due to the competition of length scales. 
Later studies have often focused on this nonuniversahty, 
e.g. in the context of wetting phenomena in supercon- 
ductors |l^, alloys ^9 20 1 and related systems pif . 

The present paper is organised as follows. In section || 
we introduce our Ginzburg-Landau free energy functional 
and obtain its wetting behaviour in the limits of infinite 
and vanishing order parameter stiffness. At intermediate 
values of the stiffness parameter the wetting behaviour is 
found flrstly via an analytical minimisation of the func- 
tional wit hin a piecewise parabolic potential approxima- 
tion (sec. II B), and then (in sec. p^I C| ) via a numerical 
minimisation of the free energy functional to obtain the 
density/order parameter proflles. In section III A we re- 
port the results of grand canonical Monte Carlo studies 
of a symmetrical binary Lennard- Jones fluid at an attrac- 
tive structureless wall. The density and order parameter 
proflles with respect to the wall are obtained along a 
sub-critical isotherm for a number of different wall-fluid 
potential strengths. Finally we compare and discuss the 



mean-field and simulation results in section [V 



II. GINZBURG-LANDAU THEORY 

Our theoretical studies are based on a generic 
Ginzburg-Landau functional for a system with two or- 
der parameters ^(r, z) and m(r, z): 



with the bulk free energy density 

f m, ffl) = -m H — -d) m H m 

J\ 2 4 2 4 

+fi(j>- Km^(j) (2.2) 
and the bare surface free energy at the wall 



drj^ rfz{|(V0)2 + |(Vm)2 + /(m, 0)} 

dffs{m,<f>)U^o (2.1) 
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(2.3) 



The z-axis is taken to be perpendicular to the waU 
and J dr integrates over the remaining spatial dimen- 
sions. In our case, the quantity m is related to the 
difference between the partial densities of the two com- 
ponents, m cx (pa — Pb), and 4> to the total density, 
(f> cx {p — po), where the reference density po is chosen 
in the liquid/vapour coexistence region such the cubic 
term proportional to (p — po)^ in (2.2) vanishes. Be- 



low the liquid/ vapour critical point, it is convenient to 
set the units of 4>, to, F and of the length such that 
bm = = = g = I, and to define 9 = am — 1- The 
bulk free energy density then takes the form 



11^ 

/(m,0) = --(/)2 + -(/,4__™2 

— + k(1 — 0)m^ . 
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(2.4) 



The bulk properties of this model have been discussed 
earlier [Q. A A-line 0\{p-) of continuous transitions sep- 
arates the mixed fluid from the dcmixed fluid at large 
negative p, corresponding to large densities 0. As long 
as K < 1, it is terminated by the onset of liquid/vapour 
coexistence at a critical end point {Ocep = 0, Peep — 0). 
The parameter p, is field like and is temperature like, 
9 <x. {T — Tcep), where T^ep is the critical end point tem- 
perature. Above 6'cep, liquid/vapour coexistence is en- 
countered at p = 0, and below 9cep, at 



Pc 



(1-k2)' 



(2.5) 



The coexisting liquid and gas phases are characterized by 
the order parameters (to linear order in p) 







in the gas phase, and 



Kp 



1 



-l-p/2 (2.6) 



= 1-| + |to;^ (2.7) 



in the liquid phase. These expressions are also valid in 
the regime where the liquid or gas phase are metastable. 

Minimizing the functional {2A) yields the Euler La- 
grange equations 



df 

ITT 7l 



^ (dz)2 ' {dzY 

with the boundary conditions 



5/ 
dm 



9 



dz 



dfs 



dm 



z=Q 



1- 



z=0 



dm 



(2.9) 



We wish to study a situation where the mixed liquid 
(to = 0) wets the wall at p = (coexistence between 
vapour and mixed liquid). To ensure this under all cir- 



cumstances, we choose = —(PoCa, with 0o > 



and 



take the limit — > oo, which is equivalent to constrain- 
ing the surface density at the fixed value 0(0) = 00- The 
surface coupling Cm is taken to be positive. It accounts 
for weakening of the dcmixing tendency at the surface 
due to the reduced number of interacting neighbors. 

One possible solution of the Euler Lagrange equations 
describes a mixed film at a wall. In this case, m{z) — 
everywhere and one is left with one order parameter 
only. The bulk value of in the metastable mixed 
phase is (p^^ = 1 — p/2. The standard way of solving the 
problem |Q] shall be sketched briefiy for future reference. 
One begins by identifying the integration constant 



1 



(2.10) 



which gives an expression for d(j)/dz as a function of (j). 
The surface free energy can then be expressed as an in- 
tegral over (j) 



(2.11) 



and the excess density (jJi'xc = \ /q°° dz \4'{z) 
the surface can be calculated via 



a(0)1 



at 



1 



As long as |/(0f -/(0_)| « (^V"^ - V"(0V" ; which 
is true for p ^ 2, the main contribution to this integral 
stems from (j) values around 0^'' . The numerator in the 
integrand can then be expanded around 

this to second order and assuming p ^ (0o 
obtains 



6(0)) 



v/2(/(0)-/(0L"^))' 



(2.12) 



(0) 

^ . Oarrymg 

^+), one 



6(0) 



ln( 



4(00 



B < cj^f) (2.13) 



•0 > 0f ). (2.14) 



Above the bulk demixing transition, the mixed film thus 



wets the wall at coexistence (p 
maintains a finite thickness for 0o < 

00 > 



0) for 00 > 0^"*, and 



We will choose 



jI^-' hereafter. From eqn. 2.1l|, one calculates the 



surface free energy to leading order in p and 



(0) ^ 2V2 . 



1 

72 



bf) 



(2.15) 



Below the bulk demixing transition, p = pc > at coex- 
istence and the thickness of the mixed film remains finite 
under all circumstances. 

In the following, we shall first analyse the wetting be- 
havior for the limiting cases where the order parameter 



2 



varies varies on very short length scales (^^jm*^ 0) 
and on long length scales (7/™+^ — > 00) compared to the 
density. Then we will discuss the general case of inter- 
mediate 7, first analytically in an approximation where 
the potential (^) is replaced by a piecewise quadratic 
potential, and then numerically with the full potential 

(il- 



A. Limiting cases 

We consider first the wetting behavior at (7/TO^^ — > 
0). In this case, m adapts locally to </>, and the order 
parameter profile miz) can be written as m(0(z)) with 
m{<^) = e + 2K{(j) - 1) for < 1 - e/2K and jn{(j)) = 
otherwise. Hence we are left with the effective one or- 
der parameter problem of calculating the density profile 
in the slightly altered potential /((/)) = f{m{(p), (f). 
Since /((/)) is a smooth function with two minima, one 
can proceed as sketched above for the mixed film, with 
the analogous result: The demixed film wets the wall at 

</)>0;. 

The analysis of the opposite case, (7/m!j_^ — s- 00), is 
somewhat more involved. Here cj) adapts locally to m\ 
however, the bulk equation df /d(f> = (p^ — ip+ii—Km'^ = 
has two solutions (j)± (m) . One conveniently separates the 
profiles into four parts (I) - (IV) as indicated in Figure 
|l|. The regions (I) and (III) are narrow slabs where (/)(z) 
varies rapidly and m can be approximated by a constant, 
m = mo at the surface (I) and 7712 at the interface (III). 
In (I), (j) drops from it's surface value (po to the local 
equilibrium value (/)_|_(mo), and in (HI), it switches from 
(/)+(m2) to (/)_(m2). The other two regions, (II) and (IV), 
are much wider; The order parameter m{z) varies slowly 
and ipi^) adjusts locally to m{z)^ such that (j) = (j)-\-{ni) 
in (II) and = 0_(m) in (IV). 

To make the argument more quantitative, we specify 
the actual subdivision of the excess free energy of (p?l|), 

J^exc = fs{mo, 0o) J^ii + Tin + Tiv (2.16) 

with 



III IV 
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FIG. 1. Schematic sketch of density and order parameter 
profiles in the limit '^jm*^ — > 00. See text for more explana- 
tion. 



Ti = dz[-{ — )^ + f{mo,(l3)-f{mo,(l3+{mo))] 

J 



^ 1 djTi 
dz[--fi — f + f{m,(b+im,))] 



'-2 " 

'^^[o(:r)^ + /("^2,0) - /(m2,0+(rn2))] 
'^^^l^)' + /("^2, 0) - /(m2, 0-(m2))] 
'rf4^7(^)' + /("i,0-(r«))]. 



The calculation for the regions (I), (HI), and (IV) can 
proceed in an analogous way as sketched earlier for the 
mixed film: The profiles of (p(z) in (I), (III), and of to(z) 
in (I V) are monotonic and the integration constant (cf. 
2.10) is known (zero). One obtains to leading order in /i 
and (00 - 0+(nT-o)) 



Fill = 
Fiv — 



2V2 

2 

TO2 



>o - 4>+{'mo)f + • 

+ 0((/i-KTO2')') 
^{Ak + K^l - 9) 



(2.17) 
(2.18) 
(2.19) 



In the region (II), the integration constant is unknown 
1 ,dm 



^7(^)^ - fim, 0+(m)) p, 



(2.20) 



with p > if the profile of m(z) is monotonic, and p < 
if m{z) is nonmonotonic, like in Fig. 0. A connection 
between p and the width I of the film can be established 
using / — J™" dm/\dm/dz\ in the first case, and 



dm 



\dm/dz\ 



dm 



\dm/ dz\ 



in the second case, where m„iax solves p = 
~ f{mmax,4'+{''f^max))- Ncxt wc expand the function 
f{m,(j)^{m)) about it's minimum m*^, leading to 



m 2 1~ 



/(m. Mm)) « (1 - ^)m+\{--^ - If - -) (2.21) 



One deduces the characteristic length scale, 
A = 



7 



1 



2(1 - k2) m* 



(2.22) 



which grows very large in the limit ^ jm*^ — cxj. The re- 
sult for Tii can therefore be expanded in powers of e~'/^. 
After adding up all contributions (I)-(IV) and minimiz- 
ing with respect to m2, the total excess free energy of 
the demixed film takes the form F^^c = Fsurfimo, (j)o) + 
Fint + V{1) with the surface contribution 



3 



Fsurf = fs{mo, (j)o) + y2^^o~ 0+("lo)) 

the interface contribution 
2V2 



8A^c 



(2.23) 



(2.24) 



and a surface/interface interaction term 

V{1) = 2(^ - ^c)l - 32A/Xc(l - mo/m;)e"'/^ 

+16A/ic(yl - (1 - mo/m;)2B)e"2//A^ (2.25) 



where B = [A^7(k - 771+2(1 - k2)/4]-i^ a,nd A = -3 or 
A = 1, depending on whether or not the profile m{z) is 
monotonic. 

The result can now be discussed. At 7779 < 771^, the 
leading term e~^^^ of the potential V{1) is attractive, 
and wetting is not possible. At 7779 > 771!^, an infinitely 
thick demixed film is metastable at coexistence. It's free 
energy difference to the mixed film AF — F^xc — F^xc is 
up to third order in 777I 



AF 



i= (VW^) ( 1 + (1 - r77o/777;)2) 777;^ 

+ (a„/V2-K((/.o-'/'+(™o)))r772) (2.26) 



00 can be taken in two ways: either 
at fixed 7. In the 



The limit 7/777*^ 
7 ^ 00 at fixed 
first case the first term in eqn. (2.26) dominates and the 
free energy of the demixed film exceeds that of the mixed 
film: The film remains mixed and dewets accordingly. 

The second case is more subtle. Here, the second 
term dominates, and the free energy of the mixed film 
may be less favorable, depending on the ratio of Cm and 
(00 — <?!)+). Note that the density enhancement at the 
surface, •\/2/t(0o ~ 0+)j acts as an additional surface cou- 
pling, which opposes the effect of Cm- The parameter 
Cm accounts for the direct reduction of interacting neigh- 
bours right at the surface. It is counterbalanced by the 
fact that the density 0o close to the surface is higher 
than in the bulk. If the latter effect dominates, the film 



demixes at the surface even for 



or T ^ T, 



cep • 



B. Analytical results in a piecewise parabolic 
potential 

At fixed 777^ , we have seen that the demixed film wets 
the substrate in the limit 7^0, where the order param- 
eter 771 varies much faster than the density </), and dewets 
at 7 ^ 00, where the density varies much faster than the 
order parameter. Now we consider intermediate values 
of 7, where the two characteristic length scales become 
comparable. Far from the critical end point, this is the 
usual case in a binary liquid, since the interaction ranges 



responsible for liquid/gas separation and demixing are 
comparable. 

In order to carry further the analytical an alys is, we 
approximate the free energy density /(</>, tr) (2.4) by a 
piecewise quadratic form 



with three pieces corresponding to the gas phase and the 
two liquid phases, separated by the lines 

(t>sep{m) = -K{m^ + - 17771777";) + [1/2. (2.28) 

and 777 = at > 0sep(O). Here /2 = /i — /ic, (J = — 1 
for (j) > (j>sep{rn) (gas phase), a = +1 for (p < (j)sep{m) 
(liquid phases), and the parabolae are adjusted to the 
leading terms in the expansion of the functional (2.4) 
about its minima. 



-1 




/ 



/2 ^ 
VO 4k J' 



(2.29) 



for (f) < (psepim) (gas phase), and 



(^) 

V 777 / 



1 + nmY/2 
m J V ±m* 



_ 2 =F2k777* 



for 



> 



(2.30) 

hsepim) (liquid phases), where the upper sign 



holds for 777 > 0, the lower for 777 < 0. The choice (2.2S) 
of (psep ensures that the potential f(m, (p) is continuous. 

In such a potential, profiles of demixed films corre- 
spond to paths in the (0, 771) space which can be sep- 
arated into three parts: (i) moving in one of the liq- 
uid regions from (0o,wo) to (0i,77i — 0); (ii) following 
the edge (771 = 0) between the two liquid regions from 
(01, 0) to (0sep(O), 0); (iii) moving in the gas region from 
(0sep(O),O) to (—1,0). On principle, a direct transition 
from (i) to (iii) is conceivable. For the parameters 0o 
of interest, however, such profiles turn out to be ener- 
getically less favorable than the profiles which have an 
intermediate (ii). Profiles of mixed films have two parts 
(ii) and (iii) only. We shall denote ^(i) = I, l^u), and ^m), 
the length of the slab spent in region (i), (ii) or (iii), 
respectively. 

At given slab length and boundary conditions, the free 
energy in each of the slabs can be calculated exactly using 



I r ,du 



2 
1 



^{(77(0) + 77(0)2 tanh — 



(2.31) 

(77(0)- 77(0)2 COthi^}. 



The calculation is straightforward in the regimes (ii) and 
(iii). In (i), the free energy functional has to be diago- 
nalized first: 



(0 











' , dw , 9 
L dz 


w 


L az 




+ 





4 



with 



1,2 



V w 




where we have defined 



2 'x-^-l' 



(2.32) 



(2.33) 



The parameter 5 or alternatively determines the 

wetting behavior. Fig. || shows the two length scales Ai 
and A2 as a function of ^/m^'^. The length Ai is always 
larger than A2. At ^jm*^ 3> 1 or 5 3> 0, it charac- 
terizes the spatial variations of m(z) and grows linearly 
with ^lra*^\ at <C 1 or (5 ^ 0, it characterizes 

the variations of <i>{^z) and remains largely independent 
of 7/™+ . These are the limiting regimes discussed in 
the previous subsection. At ^jmZ^ w 1 or (5 w 0, both 
Ai and A2 are related to linear combinations of <\){z) and 
m(z). 

The further calculation proceeds as follows: The free 
energy in (iii) is given by 

% = %/2/2(l + (0,,p(O) + l)2). (2.34) 

In the region (ii), the result for the free energy is ex- 
panded in powers of e~^'<") up to the second order and 
minimized with respect to The free energy calcu- 

lated in (i) is expanded up to second order in powers of 
g-i/^i and up to first order in e"'/'*'^, where / = The 
three contributions are then added up, and the sum is 
minimized with respect to 0i and mo at given surface 
coupling Cra- The solution has to be compared with the 
free energy of a mixed film, which is calculated analo- 
gously. 



We only report the result for the case Cm = here. 
The expressions obtained for arbitrary Cm are more com- 
plicated, but qualitatively similar. Without loss of gener- 
ality, we can assume to > in the demixed film. As long 
as Too > 0, the surface order parameter toq and the free 
energy difference Ai^ between the mixed and demixed 
film can then be expanded as 



AF(/) = AF(,,) + TO + r2e-'/^i -f T2 

Using the abbreviations 

-, and K± — 



-2l/\i 

-2//A1 
(2.35) 



A1A2 



A1A2 ' e±*'Ai -t- e=F« A2 ' 

the coefficients can be written as 

Lo = (A2-1 - A-i)X_ • (00 - 0) 
ii,2 = -2i^o^+^-e±^ • V7to; 



(2.36) 



To = Ko{K+ ■ jmf + ■ (00 - 0)') 
ri,2 = ±2i^oi^+i^-/A2,i • Vl^X ■ (00 - 0) (2.37) 
ra - KoK-Kle'{e-'/X2 - e'/AO/Aa • 7mf 

+ KaKle-\\ + 8/f+e-VAi)/A2 • (0o - 0)', 

and with h = TO!j_'*(l — k^) — 2jl, p — (j)o ~ 4> + kto^'^. 

When taking the limits S —>■ ±00, one recovers qualita- 
tively psl the behavior discussed in the previous section. 





FIG. 2. Length scales Ai and A2 vs. j/m'^^ at k = 0.5. 
Thin line shows Ai/2 for comparison. 



FIG. 3. Coefficients of interfacial potential vs. '^jra*^ 
for K = 0.5, "^o = 1-5, and surface coupling Cm = 0. (thick 
lines) and Cm = 0.2 (thin lines). 
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The function AF{1) can be conceived as an effective 
interface potential for the demixed film. The parameters 
Ti for a choice of (f>o {(tri — 1-5) and two values of Cm, 
Cm = according to (2.37) and Cm — 0.02, are shown in 



Fig. H One finds that ti is always positive, T2 is always 
negative, and changes sign from positive to negative 
as ^ jm*^ increases. The leading term of the potential 
is thus positive, and one expects a first order wetting 
transition a nd a p rewetting line. On the other hand, the 
expansions ( 2.35 ) are only valid as long as the s urfac e 
order parameter mo is positive. According to eqn. ( p. 36 ), 
the coefficients li of the expansion for mo(0 £^re negative 
except for the zeroth order term lq. Hence mp decreases 
with film thickness and may vanish at some thickness Ic- 
In this case, the film mixes continuously at lc^ and the 
prewetting line turns into a second order demixing line 
sufficiently far from coexistence. 



C. Numerical solution 
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FIG. 4. Phase diagram in the fi — <j)o plane. Parameters 
are k = 0.5, 6 = 0.1 and 7 — 1., Cm = 0.2. Solid line indicates 
first order transition, dashed line second order transition. 



The analytical results of the previous subsection pro- 
vided insight into the competition of length scales in the 
binary fiuid and the wetting scenarios which can be ex- 
pected on a wall as a result. However, a reliable calcula- 
tion of actual phase diagrams, including the details of the 
prew etting line, is not possible on the basis of the expan- 
sion ( 2.35 ). We have thus supplemented the analytical 
work by a numerical minimization of the functional (2.1) 
in the fi — (j)o plane for selected sets of parameters 7 and 

Cm- 

The problem is simplified considerably due to the fact 
that <f>{z) is a monotonic function of 2, i.e., m{z) can 
be expressed as a function m{(j)). The bulk free energy 
functional in (2.1) can thus be rewritten as 



T=l dz{i(l + 7($)^)(f )^ + /(m(<^),0)} 



2 

4>o 



d<ty ^l + 7(^)V/(m(<^),0)-/(O, </.!), 



where the integration constant ( 2.1C ) has been identified 
and exploited as usual. Minimization with respect to the 
function m{(j)) leads to the Euler-Lagrange equation 



27/(m,(/.)- 



^ ^dniy^^ ^ df dm df ^ 
uujj- d(j) dm d(j) dcf) ' 

(2.39) 

which we have solved using the Vcrlet algorithm. 

Some results are shown in Figs. ^, ^, and |^. As an- 
ticipated in the previous subsection, we find a first order 
wetting transition, a discontinuous prewetting line and a 
continuous demixing line. At surface coupling Cm > 0, 
the demixing line joins the prewetting line in a surface 
critical end point (Fig. The prewetting line separates 
a demixed thick film from a mixed thin film (see profiles 




FIG. 5. Density and order parameter profiles for the co- 
existing mixed and demixed film at the point in the /j. — (f>o 
plane indicated by the arrows in Fig. 6. 

in Fig. 1^) before reaching the critical end point, then two 
demixed films of different thickness, and finally vanishes 
in a critical point. On decreasing the surface coupling 
Cm, the critical end point and the critical point move 
closer to each other, until they merge in a surface tricrit- 
ical point. 

Fig. 1^ shows two cases of phase diagrams in the 0o ~ M 
plane for Cm — and two different 7 at fixed 9, i. e., 
at fixed bulk order parameter m^. With increasing 7, 
the prewetting line shifts towards larger (j)Q and extends 
deeper into the off-coexistence region. As 7 00, it 
moves to (/)o — + 00, the film remains mixed and thin at all 
finite 4>o. At 7 ^ 0, on the other hand, the line becomes 
fiat, approaches 4>'^, and the tricritical point where it 
turns into a second order line moves to fj,t — > fJ.c- The 
numeri cal r esults thus agree with the conclusions from 
section HA. 
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Fig. demonstrates what happens if instead of mak- 
ing 7 larger, one increases the characteristic length scale 
of order parameter fluctuations by decreasing 9, i.e., ap- 
proaching the critical end point (reducing m*^). Far from 
liquid/vapour coexistence, the transition line still moves 
towards larger (pQ. However, the effect reverses close to 
coexistence, the demixing transition is now shifted to 
lower surface densities 0o. Furthermore, the length of 
the prewetting line shrinks instead of growing. 



1.30 



1.25 



1.20 



1.15 



1.10 



1.05 ^ 
0.0000 




0.0005 
U- U 

v r CI 



0.0010 



0.0015 



FIG. 6. Phase diagrams in the /i — (jio plane for different 
7 with 9 = 0.1, surface coupling Cm = 0., and parameters as 
in Fig. 6 otherwise. Solid lines indicate first order transition, 
dashed lines second order transitions. 
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FIG. 7. Phase diagrams in the jj, ~ (j}o plane for different 
9 with 7 = 1., surface coupling Cm ~ 0. and parameters as 
in Fig. 6 otherwise. Solid lines indicate first order transition, 
dashed lines second order transitions. 



III. MONTE CARLO SIMULATIONS 

In this section we describe Monte Carlo simulation 
studies of the subcritical wetting behaviour of a sym- 
metrical binary fluid at a structureless wall. 

A. Model and simulation details 

The system we have studied is a symmetrical binary 
fluid, having interparticle interactions of the Lennard- 
Jones (LJ) form: 



/ \ 12 













(3.1) 



We made the following choice of model parameters: 

CTii = 0-22 = o'i2 = cr = 1; £11 = £22 = e; ^12 = 0.7e. 
i.e. interactions between similar species are treated iden- 
tically, but those between unlike species are weakened. 
The inter-particle potential was truncated at a distance 
of i?c = 2.5(7 and no long-range correction or potential 
shift was applied. 

The fluid was confined within a cuboidal simulation 
cell having dimensions x Py x D, in the x,y and z 
coordinate directions respectively, with P^ — Py = P. 
The simulation cell was divided into cubic sub-cells (of 
size the cutoff Rc) in order to aid identification of particle 
interactions. Thus P = pRc and D = dRc, with p and d 
both integers. To approximate a semi-infinite geometry, 
periodic boundary conditions were applied in the x and y 
directions, while hard walls were applied in the unique z 
direction at z = and z — D. The hard wall at 2 = was 
made attractive, using a potential designed to mimic the 
long-ranged dispersion forces between the wall and the 
fluid pH: 



Viz) 



(3.2) 



z measures the perpendicular distance from the 
is a 'well-depth' controlling the interaction 



Here 
wall, 

strength, and we set cjyj — 1. No cutoff was employed 
and the wall potential was made to act equally on both 
particle species. 

Monte-Carlo simulations of this system were performed 
using a Metropolis algorithm within the grand canoni- 
cal (^,y, r) ensemble |25). Three types of Monte-Carlo 
moves were employed: 

1. Particle displacements 

2. Particle insertions and deletions 

3. Particle identity swaps: 1 — > 2 or 2 1 

To maintain the symmetry of the model, the chemical po- 
tentials /ii and /i2 of the two components were set equal 
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at all times. Thus only one free parameter, = /ii = /i2, 
couples to the overall number density p = {Ni + N2)/V. 
The other variables used to explore the wetting phase 
diagram were the reduced well depth e/ksT and the re- 
duced wall potential e^/kBT. During the simulations, 
the observables monitored were the total particle density 
profile 

p(z) = [iVl(z) + iV2(z)]/P^ (3.3) 

the number difference order parameter profile, 

n{z)^\Ni{z)-N2{z)\/P'' (3.4) 

These profiles was accumulated in the form of a his- 
togram. Other observables monitored were the total in- 
terparticle energy and the wall interaction energy. 

The choice of system size was, as ever, a compro- 
mise between minimising finite-size effects and maximis- 
ing computational throughput. Tests showed the profiles 
to be largely insensitive to the size of the wall area and 
hence p = 7 was used, this being the largest computation- 
ally tractable size consistent with the necessary choice of 
the slit width d. The latter must clearly be considerably 
larger than the film thicknesses of interest in order to 
prevent the liquid film directly interacting with the hard 
wall at z = D. In the results presented below, the typ- 
ical slit width used was d — 16, corresponding to some 
40 molecular diameters. For thin films a narrower slit of 
width c? = 8 was used. 



B. Wetting behaviour along a subcritical isotherm 

Accurate knowledge of bulk phase behaviour is 
an essential prerequisite for detailed studies of near- 
coexistence wetting properties. In the present model, the 
locus of the liquid vapour coexistence curve and location 
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FIG. 8. The phase diagram in the ji-T plane of the sym- 
metrical binary Lennard- Jones fluid model described in the 
text. Also shown is the location of the critical end point and 
the isotherm along which the wetting properties were studied 



of the critical end point are already known to high pre- 
cision from a previous MC simulation study The 
phase diagram in theu-T plane (in standard Lennard- 
Jones reduced units [^) is shown in fig. ||. The critical 
end point is located at Tcep — 0.958(3), ^cep — —3.017(3) 
[^,|| . We note that although the locus of the coexistence 
curve is known to five significant figures, the position of 
the CEP along this tightly determined line is known only 
to three significant figures. 

To determine the wetting properties at temperatures 
below Tcep, the number density profile p{z) was stud- 
ied along the isotherm T = 0.9467 as coexistence was 
approached from the vapour side. To achieve this, the 
chemical potential was incremented up to its coexistence 
value Pcx{T) in a sequence of 6-10 steps of constant size 
A/x = 0.0025. This procedure was repeated for a num- 
ber of different values of the wall-fluid potential strength 
twi allowing the influence of this parameter on the wet- 
ting behaviour to be ascertained. In all, six values of 
the e„ were studied (e„ = 1.0,1.7,1.75,2.0,3.0,4.0). We 
describe the wetting behaviour for each in turn. 

For e^j — 1.0, fig. ^ shows that although the film thick- 
ness grows very slightly as coexistence is approached, it 
never exceed two molecular diameters. At no point in the 
profile does the density attain that of the liquid phase 
(p « 0.6). The presence of a thin wetting layer right up 
to coexistence implies incomplete (partial) wetting. 

Increasing the wall potential to — 1-70 [fig. 0, 
results in considerably more structure near the wall com- 
pared to tw = 1-0, with clear density oscillations arising 
from excluded- volume 'packing effects' The profile 

is much more responsive to changes in the chemical po- 
tential and reaches a thickness of 4-5 molecular diameters 
close to coexistence. 

For e^, = 1.75, however, the situation changes quali- 
tatively, as shown in fig. 1 1 , On increasing the chemical 
potential, a clear jump is observed in both the thickness 
of the film, and the value of its density. In the thick film, 
the density of a significant portion of the film is that of 




FIG. 9. Density profiles for ~ 1.0. Data are shown for 
7 values of /i — fj,cx in the range [0,-1.5] 
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FIG. 10. Density profiles for = 1.7. Data are shown 
for 8 values of /i — ficx in the range [0, —1.525] 
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FIG. 11. Density profiles for e„ = 1.75. Data are shown 
for 6 values of /x — Hcx in the range [—0.025, —1.5] 
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FIG. 12. Density profiles for e„ = 2.0. Data are shown 
for 8 values of /i — jj,cx in the range [—0.025, —1.6] 



the bulk liquid. This thin-thick jump constitutes a 
prewetting transition, as previously observed in simula- 
tion studies of lattice gas models , Lennard- Jones flu- 
ids ||l^,|lj,^,^ as well as experimentally [^ . 

As the wall potential is increased to — 2.0 (fig. p^ , 
the sharp prewetting transition is lost and instead the 
film thickness increases smoothly as [i approaches its co- 
existence value. This suggests that here the system is 
above the prewetting critical point ( ]2^). 

On increasing to 3.0, a new feature emerges 
(fig. [l3|(a)). As the chemical potential increases, the 
thickness of the film initially increases smoothly with in- 
creasing [I. However, once the thickness reaches some 10 
molecular diameters, a large jump occurs to a thickness 
of about 15 molecular diameters. Concomitant with this 
jump is a demixing of the film as a whole, as seen in the 
order parameter profile fig. [l^(b). The size of the jump in 
the layer thickness appears to reduce as the wall strength 
is increased to e^, = 4.0 (fig. suggesting a weakening 
of the transition. 
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FIG. 13. (a) Density profiles for = 3.0. Data are shown 
for 6 values of p — \iax in the range [—0.025, —1.55] 

(b) The corresponding order parameter profiles n(z). 
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FIG. 14. Density profiles for — 4.0. Data are shown 
for 8 values of /i — jj,cx in the range [—0.025, —1.6] 

IV. DISCUSSION 

The Monte Carlo simulation results at subcritical tem- 
peratures provide evidence that the mean field calcula- 
tions correctly identify the qualitative wetting behaviour. 
They show that depending on the fluid-wall interaction 
strength e^,, a number of different wetting scenarios oc- 
cur as liquid-vapour coexistence is approached from the 
vapour side. At small e^, only a very thin film builds 
up on the wall. For intermediate values of e^;, a first 
prewetting transition is observed from a thin mixed film 
to a thick liquidlike mixed layer. Further increasing 
induces a second prewetting transition between a mixed 
liquidlike layer and a thicker demixed filni, the situation 
being very similar to that shown in figure |5|. The abrupt, 
first order, character of this latter transition appears to 
weaken on further increasing in accord with the the- 
oretical predictions. 

We will now attempt to set our results within the con- 
text of the bulk phase diagram of the binary liquid. To 
this end, we discuss the possible wetting scenarios in the 
vicinity of the critical end point Tcep- As previously ar- 
gued in the introduction, for temperatures T < Tcep suf- 
ficiently close to Tcep, the bulk correlation length ^ of the 
demixed liquid is larger than the thickness l^, of a mixed 
liquid layer at the wall. The state of order of the film 
thus depends strongly on the boundary conditions of the 
two interfaces confining the liquid layer. The nonselec- 
tive liquid-vapor interface always favors mixing due to 
the reduced number of interacting neighbors in the inter- 
facial region. The liquid-substrate interface, on the other 
hand, can either favour mixing or demixing depending on 
the strength of the fluid- wall potential. For a weakly at- 
tractive wall potential, mixing is favoured because the 
particle density at the wall is low and the presence of the 
wall reduces the number of interacting neighbours. For a 
strongly attractive wall, however, the high density at the 
wall can counteract the missing neighbour effect leading 
to an overall demixing tendency. 




Density p 




Density p 




Density p 



FIG. 15. Some possible schematic wetting phase diagrams 
in the temperature-density plane, (a) Weakly attractive wall: 
Critical end point Tcep of the A-line is a critical wetting point, 
below which the wall is not wetted by the liquid, (b) Interme- 
diate attraction: Demixing induced first order wetting transi- 
tion at T < Tcep with prewetting line which evolves into a sec- 
ond order demixing line, (c) Strong attraction: Complete wet- 
ting at coexistence everywhere, but detached prewetting line 
or continuous demixing transitions ofi' coexistence (in films 
of finite thickness). Hatched area indicates the possibility of 
conventional wetting transitions at lower temperatures. 
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If the net effect favours mixing at the wall, a contin- 
uous demixing of the layer as coexistence is approached 
can be excluded. A first order transition involving a dis- 
continuous increase of the film thickness upon demixi ng is 
still conceivable. However, we have shown in section [I A, 



that (at the mean field level, at least) the demixed wet- 
ting film has a higher free energy than the corresponding 
mixed film provided the correlation length of composition 
fluctuations is sufficiently large. 

At walls which suppress demixing, the film is thus al- 
ways mixed close to the critical end point, and its thick- 
ness below the critical end point is finite. Hence the 
critical end point is automatically a critical wetting point. 
The resulting phase diagram is shown schematically in 
fig. |l^ (a). Note that the wetting transition here is 
pinned by a bulk phase transition, a situation somewhat 
reminiscent of triple-point wetting |3^Jll[] . 

The situation changes if the substrate favours demix- 
ing. In this situation, one component segregates to the 
surface of the film already slightly above Tcep, and the 
order propagates continuously into the bulk of the film 
at Tcep- The fi lm remains wet at T^ep- From the results 



of section II C (in particular. Fig. 0), one can deduce two 
possible scenarios. The film may still exhibit a first or- 
der wetting transition to a nonwet state at a temperature 
below Tcep (e.g. in Fig. 0) at (f>o — 1.14). The discontin- 
uous phase transition at liquid/ vapour coexistence then 
spawns a prewetting line which eventually switches into 
a second order demixing line and loops around the crit- 
ical end point as suggested in Fig. |lj (b). If the wall 
is strongly attractive (e.g. at (j)o = 1.27) in Fig. |^), 
the wall wets at all temperatures, the prewetting line 
detaches from the coexistence line and is continued by 
second order demixing lines both at the high and low 
temperature side as sketched in Fig. |l^ (c). 



Acknowledgements 

NEW thanks the Royal Society of Edinburgh, the EP- 
SRC (grant no. GR/L91412) and the British Council for 
financial support. 



[1] L.D. Gelb, K.E. Gubbins, R. Radhakrishnan, M. 
Sliwinska-Bartkowiak, Rep. Prog. Phys. 62, 1573 (1999). 

[2] For reviews on wetting see P. G. de Gennes, Rev. Mod. 
Phys. 57, 827 (1985); M. Schick, in Les Houches, Ses- 
sion XLVIII - Liquids at Interfaces, J. Charvolin, J. F. 
Joanny, and J. Zinn- Justin eds. (Elsevier Science Pub- 
lishers B.V., 1990); and Ref. 

[3] M. N. Barber, in Phase Transitions and Critical Phe- 
nomena, Vol. 8, p. 145, C. Domb, J. L. Lebowitz eds. 
(Academic Press, New York, 1983); V. Privman in Fi- 



nite Size Scaling and Numerical Simulation of Statistical 
Systems, V. Privman ed. (World Scientific, Singapore, 
1990). 

[4] M.P. Nightingale and J.O. Indekeu, Phys. Rev. Lett. 54, 

1824 (1985); ibid 55, 1700 (1985); M. Krech, Phys. Rev. 

E 56, 1642 (1997); M. Krech, J. Phys.: Condens. Matter 

11, R391 (1999). 
[5] N. B. Wilding, M. Krech, Phys. Rev. E 57, 5795 (1998). 
[6] A. Mukhopadhyay and B.M. Law, Phys. Rev. Lett. 83, 

772 (1999). 

[7] N. B. Wilding, F. Schmid, P. Nielaba, Phys. Rev. E 58, 
2201 (1998). 

[8] M. E. Fisher, P. J. Upton, Phys. Rev. Lett. 65, 2402 
(1990); M. E. Fisher, M. C. Barbosa, Phys. Rev. B 43, 
11177 (1991). 

[9] N. B. Wilding, Phys. Rev. Lett. 78, 1488 (1997); N. B. 

Wilding, Phys. Rev. E55, E55, 6624 (1997). 
[10] S. Dietrich in Phase Transitions and Critical Phenomena, 

Vol. 12, p. 1, C. Domb and J. L. Lebowitz eds (Academic 

Press, New York, 1988). 
[11] S. Dietrich, M. Schick, Phys. Rev. B 33, 4952 (1985). 
[12] I. Hadjiagapiou, R. Evans, Mol. Phys. 54, 383 (1985). 
[13] E. Kierlik, M. L. Rosinberg, Y. Fan, P. Monson, J. Chem. 

Phys. 101, 10947 (1994). 
[14] Y. Fan, J.E. Finn and P.A. Monson, J. Chem. Phys. 99, 

8238 (1993). 

[15] A. Plech, U. Klemradt, M. Huber, J. Peisl, Europhys. 
Lett., in press. 

[16] See e.g. M.J. P. Nijmeijer and J.J. Weis in Annual Reviews 
of Computational Physics IV, World Scientific, Singapore 

(1996) , and references therein. 

[17] E.H. Hauge, Phys. Rev. B 33, 3322 (1986). 

[18] J. M. J. Leeuwen, E. H. Hauge, J. Stat. Phys. 87, 1335 

(1997) ; F. Clarysse, J. O. Indekeu, Physica A 251, 70 

(1998) . 

[19] D. M. KroU, G. Gompper, Phys. Rev. B 36, 7078 (1987); 
G. Gompper, D. M. Kroll, Phys. Rev. B 38, 459 (1988). 

[20] F. F. Haas, F. Schmid, K. Binder, Phys. Rev. B, in print. 

[21] C. J. Walden, B. L. Gyorffy, A. O. Parry, Phys. Rev. B 
42, 798 (1990). 

[22] S. Dietrich, M. Schick, Surf. Sci. 382, 178 (1997). 

[23] Quantitative difference, e.g., in the prefactors, are due 
to the fact that the underlying potentials f{m, (f>) differ 
from each other. 

[24] J. Israelachvili. Intermolecular and surface forces (Aca- 
demic Press, London, 1992). 

[25] D. Frenkel and B. Smit, Understanding Molecular Simu- 
lation Academic Press, London (1996). 

[26] D. Nicolaides and R. Evans, Phys. Rev. Lett. 63, 778 
(1989). 

[27] J.E. Finn and P.A. Monson, Phys. Rev. A39, 6402 
(1989). 

[28] Y. Fan and P.A. Monson, J. Chem. Phys. 99, 6897 
(1993). 

[29] H. Kellay, D. Bonn, J. Meunier, Phys. Rev. Lett. 71 2607, 
(1993). 

[30] R. Pandit, M. E. Fisher, Phys. Rev. Lett. 51, 1772 
(1983); J. Krim, J. G. Dash, J. Suzanne, Phys. Rev. Lett. 
52, 640 (1984). 



11 



